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ABSTRACT 

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter 
halos, two each at the mass of classical dwarf galaxies (Mvir — 10^° Mq) and ultra¬ 
faint galaxies (Mvir — 10® Mq), and with two feedback implementations. The resulting 
central galaxies lie on an extrapolated abundance matching relation from M* ~ 10® to 
10"^ Mq without a break. Every host is filled with subhalos, many of which form stars. 
Our dwarfs with M* ~ 10® Mq each have 1 — 2 well-resolved satellites with M* = 
3 — 200 X 10^ Mq. Even our isolated ultra-faint galaxies have star-forming subhalos. If 
this is representative, dwarf galaxies throughout the universe should commonly host 
tiny satellite galaxies of their own. We combine our results with the ELVIS simulations 
to show that targeting ~ 50 kpc regions around nearby isolated dwarfs could increase 
the chances of discovering ultra-faint galaxies by ~ 35% compared to random halo 
pointings, and specifically identify the region around the Phoenix dwarf galaxy as a 
good potential target. 

The well-resolved ultra-faint galaxies in our simulations (M* ~ 3 — 30 x 10® Mq) 
form within Mpeak — 0.5 — 3 x 10® Mq halos. Each has a uniformly ancient stellar 
population (> 10 Gyr) owing to reionization-related quenching. More massive systems, 
in contrast, all have late-time star formation. Our results suggest that M^aio — 5 x 
10® Mq is a probable dividing line between halos hosting reionization “fossils” and 
those hosting dwarfs that can continue to form stars in isolation after reionization. 

Key words: galaxies: dwarf - galaxies: star formation 


1 INTRODUCTION 


If ACDM is correct, then all dark matter halos hosting galax¬ 
ies - from those hosting dwarfs to those hosting giant clus¬ 
ters - should be filled with substructure ([Moore et al.|1999| 


Klypin et al. 1999 Zentner & Bullock 20031. Dark-matter 


only simulations over a vast range of particle masses and 
physical scales show that substructure persists down to the 
resolution limit of ACDM simulations ([Diemand et al.|2008 


Madau et al.|200^|Springel et aI.|2008||Boylan-Kolchin et al.| 

2009| [Klypin, Trujillo-Gomez fc Primack 20111. The ba¬ 
sic expectation is that the mass function of subhalos rises 
steadily to masses well below the molecular cooling limit of 


Mhaio ~ 10® Mq ( Tegmark et al.||l997 l, with thousands of 
sites for potentially star-forming satellites. 

Observations of our own Milky Way, on the other hand, 
have revealed the presence of only ~ 30 confirmed satel¬ 
lite galaxies (Willman et al. 2011[ McConnachie|[20T2 Be¬ 
lokurov et al. 20101, the faintest of which were all discov- 


ered by the Sloan Digital Sky Survey (SDSS) ( 

Willman et al. 

2005a|b| Grillmair|2006| Sakamoto & Hasegawa|2006| 

Irwin 

et al. |2007 Walsh, Jerjen & Willman||20071 |Grillmair 

2009 

Zucker et al. 2006a|bl Belokurov et al. 20061 |2007| 

2008 

2009 20101. However, the SDSS covers only a fraction of 
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the sky and is incomplete to the most distant and faintest 
satellites of the Milky Way. Any solution to the mismatch 
between the predicted abundance of subhalos and the ob¬ 
served counts of satellite galaxies around our own Milky 
Way will likely involve the discovery of ultra-faint satel- 


© 2015 RAS 




































































2 C. Wheeler et al. 


lites in unprobed regions of the sky, at large distances, and 
with surface brightnesses low enough that they lie just out- 


spheroidals (dSphs) show that, despite having luminosities 
that vary over nearly five orders of magnitude, they have 


side current detection limits (Benson et al. 2002a 

Ricotti] almost the same central densities, comparable to dark mat- 

& Gnedin||2005| [Koposov et al.||2008| [Tollerud et a 

.[[2008[ ter halos with Mj 

leak ~ 3 X 10® Mq (Strigari et al. 

2008 

Bullock et al.|2010 

1 . Over the last several years, searches by Geha et al. 2009 

Wolf et al.[|2010l. This is surprising, not 


PanSTARRS and VST ATLAS have failed to lend critical 
support to this idea - finding far fewer dwarf galaxy satel¬ 
lites than expected (Laevens et al.|2014 Martin et al.|2014 1. 


However, the recent discovery of up to nine new ultra-faint 
satellites in the southern sky by the Dark Energy Survey 


only because this is well above the mass where physical 
processes are believed to prevent star formation, but be¬ 
cause the mass function of potentially star-forming subha¬ 
los should be heavily populated by smaller halos. Segue 2 
may be the first of these long-searched-for small subhalos 


(DES, The DES Collaboration et al. 2015 Koposov et al. 
|2015[), the detection of three more faint dwarf satellite can¬ 


to be identified (Kirby et al. 20131, as it appears to have 
a total dark matter mass < 10“ Mq. However, based on its 


didates (Laevens et al.|[2015 Martin et al.|20f5 Kim et al. metallicity, Kirby et al. (20131 hypothesize that Segue 2 is a 


20151, and ongoing efforts to discover dwarfs at large dis¬ 
tances ( [Tollerud et al.||2015| ) provide an exciting glimpse 
into the near future and point to a much larger population 
of as-yet-undiscovered dwarf galaxies in the Local Volume. 

Interestingly, most of the newest Milky Way satellite 
candidates were all discovered near the Large Magellanic 
Cloud (LMC). While this could be nothing more than a 
location-based selection effect (the LMC happens to lie next 
to the region probed by the year one data release of DES), 
|Deason et al.| ( |2015[ ) show that several of the candidates 
could very well be satellites of the LMC, and their discovery 
highlights the potential for discovering ultra-faint satellites 
of other more massive dwarf satellites, or of isolated dwarf 
galaxies in the Local Group. The scale-free nature of the 
subhalo mass function in ACDM suggests that groups of 
subhalos should be common ([Moore et al.|1999||Li fc Helmij 


2008 Wetzel, Deason fc Garrison-Kimmel[ 20151. Because 


low-mass halos form earlier, are denser, and fall into smaller 
hosts before larger ones ( [Navarro, Frenk fc White]|1997[ ), 
it is likely that satellites of satellites or of lower mass iso¬ 
lated halos may have survived longer than their counterparts 
that fell directly into the Milky Way ( jDiemand et al.|2008| ). 
This suggests that one way to search for ultra-faint galaxies 
might be as satellites of known dwarf galaxies. Associations 
and pairs of satellite galaxies in the Local Group have been 


observed for some time ( 

Tully et al.|2006 D’Onghia & Lake 

2008 Ibata et al. 2013 

Fattahi et al. 20131, and several 


Milky Way satellites are suspected to host their own satel¬ 
lites ( Pace et al.|2014 Deason et al.|20l4 |. 

Once a group of satellites falls into the Milky Way, 
however, tidal forces will eventually disassemble the group 


and wipe out evidence of coherent structure (Sales et al. 
|2007| [Deason et al.|2015| |. Isolated massive dwarfs in the Lo¬ 
cal Group may serve as complementary targets in the hunt 
for ever-fainter dwarfs (which presumably probe the lowest 
mass dark matter subhalos). [Bovill fc Ricottij ( |2011b| ) sug¬ 
gest that even already-merged satellites of isolated dwarfs 
could be detected as “ghost halos” of ancient stellar popula¬ 
tions surrounding their hosts. Despite the obvious challenge 
of detecting ultra-faints at large distances, isolated dwarfs 
also make efficient use of telescope pointings, removing some 
of the chance inherent in any random ultra-faint search with¬ 
out losing any of the search volume. In the near future, large- 
area, deep-sky surveys such as LSST, DES, PanSTARRS, 
and SkyMapper will push current detection limits to lower 
surface brightness and may be able to see ultra-faint satel¬ 


lites of isolated dwarf galaxies if they exist (Ivezic et al. 2008 
[Kaiser et al.||200^ [Keller et al.|2007| ). 

Dynamical measurements of most Milky Way dwarf 


bare remnant of a much larger galaxy that has been severely 
tidally stripped. The common central density for observed 
dSphs might signify a low-mass cutoff in galaxy formation. 
Halos with slightly lower central densities may have been 
unable to shield themselves from the reionizing background, 
perhaps due to having a lower baryon fractions than their 
denser counterparts ( Milosavljevic fc Bromm]|2014 1 . Alter¬ 
natively, [Bullock et al.[ ( [2010[ hereafter BIO) suggest that 
the apparent common mass scale is a selection effect, and 
that a large population of unobserved “stealth galaxies” may 
reside in halos with masses just below those that host ob¬ 
served dwarfs. These low-mass halos have shallow potential 
wells, and therefore the galaxies that form within them have 
larger effective radii and lower surface brightnesses, allowing 


them to more easily avoid detection (Kaufmann, Wheeler & 
Bullock|2007 Bullock et al.|20l0 Bovill fc Ricotti|2011a|b 1. 

The number of these stealth galaxies expected to exist 
is sensitive to the presence of a low-mass cutoff in galaxy 
formation, as the lowest mass halos are expected to host 
the “puffiest” galaxies in this picture. The heating of ac¬ 
creted gas or the prevention of gas accretion by the ambient 
ionizing UV background - created when the first galaxies 
formed - can arrest or prevent star formation in the smallest 
dark matter halos, suggesting a threshold in halo mass be- 


low which all halos remain dark (Efstathiou 1992 Somerville 

2002[ Benson et al. 

2002b[ [Kravtsov, Gnedin & Klypin|2004[ 

Ricotti & Gnedin 

2005| [Moore et al.[[2006| [Strigari et al.| 

2007 

Simon & Geha[[2007 

Madau et al.[[2008 

1. It is possi- 


ble, however, that some very low mass halos were able to 


form some Population H stars before reionization (Bullock, 
Kravtsov fc Weinberg[[?000 1. If star formation ceased after 
this time, it would cause them to have uniformly ancient 
stellar populations ( Ricotti fc Gnedin[2005 1 with extremely 
low metallicity (Bovill & Ricotti 2009 Ricotti 20101. Re¬ 


markably, this appears to be the case with the known ultra 
faint dwarfs of the Milky Way ( Brown et al.[[2014 1. 

Several authors have shown that photoheating prevents 
gas from condensing in halos with masses that lie just be¬ 
low or at the common halo mass scale for Milky Way satel- 


lites ( 

Okamoto & Frenk|2009 

Nickerson et al.|2011 

Sawala 

et al. 

2014 Shen et ah] 2014 

1 . However, the specific tim- 


ing of the onset and end of reionization, the spectrum of 
the background radiation, the ability of gas to temporarily 
shield itself from these high energy photons, and particularly 
the mass resolution, can have a large effect on the number 


and minimum halo mass of galaxies in any simulation (Efs- 


tathiou|1992| [Dijkstra et al.|2004||Hoeft et al.|2006||0noi^ 


et al.|2015 1. Moreover, prescriptions for how stars and super¬ 
novae return energy back to the interstellar medium (ISM) 
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Isolated ultrafaints and tiny satellites 3 


will affect the formation history of small galaxies in the pres¬ 
ence of an ionizing background. 

In what follows, we present a series of simulations of 
galaxy formation within small dark matter halos run with 
the Gizmo code ( Hopkin^|2014 \ in “PSPH-mode” (Hopkins 
20131. Gizmo implements the Feedback in Realistic Envi¬ 


ronments (FIRE) (Hopkins et al. 20141 feedback scheme 
for converting gas into stars and capturing the energy fed 
back from those stars into the surrounding medium. Using 
these simulations, we find that the subhalos of the halos 
that surround M* ~ 10® Mq dwarf galaxies should form 
stars fairly abundantly, and produce potentially observable 
ultra-faint satellite galaxies of known classical dwarfs in the 
Local Group. Further, the ultra-faint dwarfs in our simula¬ 
tions with M* = 3 — 30 X 10® Mq that form in isolation or 
as satellites within subhalos of Mpeak — 0.5 — 3 x 10® Mq 
all have completely ancient stellar populations, as is seen 
for the known population of the Milky Way. They also have 
low-surface brightnesses, as expected for “stealth galaxies” 
forming in halos of these low masses, but may be within the 
detection capabilities of future surveys. 

This paper is organized as follows: in Sectionwe give 
the details of our simulations, including halo finding and 
merger-tree analysis. Section [^outlines our main results, in¬ 
cluding the M*—Mhaio relation for our galaxies and their star 
formation histories. In Section]^ we use the Exploring the 
Local Volume in Simulations (ELVIS) suite of dark-matter- 
only Local Group simulations ( [Garrison-Kimmel et al.|2014| 
to determine how common these objects should be and the 
likelihood of their detection. We compare to several recent 
works in Section and summarize our results and conclude 
in Section]^ 


2 SIMULATIONS 

We present six cosmological zoom-in simulations of four iso¬ 
lated dwarf galaxy halos using two implementations of sub¬ 
grid physics. The four host halos are Dwarf 1 and Dwarf 2, 
which have 2 = 0 virial masses of Mvir — 10®® Mq; and 
UFD 1 and UFD 2, which have Mvir — 10® Mq (see Ta¬ 
ble 1 for precise numbers). The Dwarf 2 halo was simu¬ 
lated three times - two runs varied the specific subgrid 
feedback implementation, while the third used a different 
choice of gravitational softening (see below), resulting in 
changes in the star formation histories of the main galax¬ 
ies (see [Onorbe et 5)]|2015| for details). We refer to these 
three cases as Dwarf 2Eariy (run with the fiducial parame¬ 
ters - the same parameters as Dwarf 1, UFD 1, and UFD 2), 
Dwarf 2Middie, and Dwarf 2Late. The subscripts refer to when 
the halos formed their stars. Every run has a dark matter 
force softening of 25 to 35 pc and a dark matter particle 
mass of m®™ = 1.26 x 10® Mq except UFD 2, which uses 
mp™ = 2.46 X 10® Mq. The UFD 1 and three Dwarf 2 sim¬ 
ulations are identical to th e “Ul trafaint” and Dwarf runs 

i2oi^.n 


Onorbe et al. 


presented by 

The initial conditions were generated using MUSIC 
(Hahn fc Abel|[2011 1, and selected from a cosmological box 


run at low resolution to 2 = 0. Dwarf 1, Dwarf 2, and UFD 1 
were all selected from 5 /i“® Mpc boxes to have typical values 
of spin parameter A, concentration, and formation time for 
their mass range, an d also to have small Lagrangian volumes 
{Onorbe et al. 20141. UFD 2 was selected from a 25 Mpc 


box and required to have no other halos of 50% or more of 
its mass within 4 f?vir at 2 = 0 and a small Lagrangian 
volume. All dwarfs in this work are isolated. Once the ha¬ 
los are identified, the particles are traced back in time and 
an enclosing Lagrangian volume is chosen and re-simulated 
at higher resolution with dark matter and gas, buffered by 
dark matter-only regions of increasing particle mass. This 
process is done according to the zoom-in techniques out¬ 


lined by Katz & White (19931; Onorbe et al. (20141 and 


with the goal of minimizing low resolution particles within 
the halos of interest. All simulations begin at a redshift of 
2 = 125 and have 0% contamination from low resolution 
particles within a distance of 1.6 Rvir at all redshifts. 

All simulations were run using the fully conservative 


cosmological hydrodynamic code Gizmo (Hopkins 20141 in 


“PSPH-mode”. This code adopts the Lagrangian “pressure- 
entropy” formulation of the smooth particle hydrodynamic 
(SPH) equations, ameliorating previous difficulties SPH 
codes had in modeling multiphase fluids ( [Agertz et al.|20(i7| . 
For feedback and star formation we use the Feedback 


in Realistic Environments (FIRE) scheme (Hopkins et al. 


20141. FIRE tracks momentum imparted locally from stel¬ 


lar radiation pressure, radiation pressure on larger scales via 
the light that escapes star-forming regions, HII photoion¬ 
ization heating, supernovae (Type I and H) heating, mo¬ 
mentum and mass loss and stellar winds from 0-type and 
AGB stars. Each star particle has an age determined by its 
formation time, an initial mass function (IMF) taken from 


Kroupa (2002 \, and a metallicity inherited from its parent 


^ Note that Dwarf 2Earl v and UFD 1 were also previously pre¬ 
sented in Hopkins et al. I 2014 1 as m9 and mlO respectively. 


gas particle. The star particle then loses mass and creates 
metals according to the STARBURST99 stellar population syn¬ 
thesis model ( Leitherer et al.||1999 l and the input IMF. Our 
simulations use an ambient ionizing UV background from 
|Faucher-Giguere et al.| ( f2009[ ), which starts at 2 = 10.65 and 
completes reionization by 2 ~ 6. 

The FIRE implementation has as few free parameters 
as possible, with star formation in the multiphase ISM nat¬ 
urally self-regulating. Giant molecular clouds (GMCs) form, 
and then are subsequently heated and disrupted by the stars 
that form within the self-gravitating, molecular gas. In order 
to ensure that the gas within the GMGs reaches a density 
high enough to form stars (n > 100 cm“®), the code requires 
extremely high mass and spatial resolution. Every run uses 
a gas particle mass of mf,^ = 255 Mq except for UFD 2, 
which uses = 499 Mq. The gas force resolution varies 
from = 1.0 — 2.8 pc in all runs except one version of 
Dwarf 2 (see below). 

The Dwarf 2Eariy run uses the fiducial subgrid param¬ 
eters. In Dwarf 2Middie and Dwarf 2Late, mass, momentum, 
and energy are deposited to particles within the SPH kernel 
according to a mass-weighting scheme, rather than volume¬ 
weighting as in the fiducial runs. In addition to this change. 
Dwarf 2Middie was run to test equal softening lengths for the 
dark matter and the gas (25 pc). All three runs of Dwarf 2 
end up with stellar masses in the central galaxies that dif¬ 
fer by only ~ 25%, suggesting that the major differences 
between the runs are more a result of stochasticity than 
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4 C. Wheeler et al. 



Figure 1. Projected (x — y plane) visualizations for all of the runs presented in this work. From left to right, the top three panels 
visualize Dwarf 1, UFD 1 and UFD 2, while Dwarf 2Barly i Dwarf 2iv[iddle s-nd Dwarf 2Late are shown in the lower panels. The dark matter 
distribution is shown in greyscale, while the stellar density is pictured in color for all subhalos with M* > 3 X 10® Mq (in at least one 
run in the case of Dwarf 2; see text for details). M* is calculated as the stellar mass within the inner 3 (1.5) kpc of the central (satellite) 
galaxy. The white line in the lower left of each panel represents 5 kpc. 


the changes to the feedback implementation and softening 
length (see Section 2.11. 


We use the Amiga Halo Finder (AHF, Knollmann & 
Knebe|2009 1 to identify gravitationally bound dark matter 


stellar and gas particles In each snapshot of the simulation. 
Because pure particle-matching algorithms generally under¬ 
perform relative to more sophisticated methods in terms of 
detection and removal of spurious interlopers between snap¬ 
shots ( Srisawat et al.|[2013 l, we track particle overdensities 
between simulation snapshots using the consistent-trees 
software | Behroozi et al.|2013b I. This is particularly impor¬ 
tant given our focus on subhalos. Because the current sta¬ 
ble version of consistent-trees - which makes use of the 
RQCKSTAR (Behroozi, Wechsler & Wu 20131 halo finder - does 
not track baryonic particles, we build a pipeline between AHF 
and ROCKSTAR that allows for a seamless transition between 
the AHF halo catalogs and the merger trees. This pipeline 
uses the halo phase space information combined with the 
dark matter particle IDs to match each halo in the cata¬ 
logue to a consistent-trees merger tree. ROCKSTAR and 
consistent-trees were used to determine the peak virial 
mass and infall times for subhalos. All other halo, galaxy 
and stellar quantities were obtained from the AHF halo cat¬ 
alogs or the raw particle files. 

All runs were initialized and simulated assuming the 
WMAP-7 cosmology as = 0.801, Ha = 0.734, flm = 0.266, 


S2(, = 0.0449, Us = 0.963 and h = 0.71 (Komatsu et al. 


20111. Further detail on the FIRE/Gizmo code can be found 
Hopkins et al. (20141; Hopkins] (20141. 


2.1 Sample Details 

Figure]^ shows the dark matter distributions of all six runs 
in greyscale with the stellar mass density overlaid in color. 
The physical scales of the panels are identical: the white line 
in the lower left corner of each visualization indicates 5 kpc. 

For our detailed analysis, we focus on objects with 
M* > 3 X 10® Mq in order to avoid effects from spuri¬ 
ous star formation at our resolution limit. The gas density 
required to form stars in our simulations is quite high and 
can only occur in dense, molecular, self-gravitating regions. 
This onerous minimal criteria suggests that the formation of 
this minimum stellar mass represents a physical star forma¬ 
tion episode. Furthermore, we find that below a subhalo 


® All galaxies that meet our minimum M* requirement have at 
least 16 star particles and had > 3 X 10^ dark matter particles 
at peak mass. All but one had > 2 X 10® dark matter particles 
at peak mass, corresponding to typical initial baryonic particle 
numbers of ~ 3 — 4.5 X 10® before feedback dispelled most of 
these. 
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Dwarf 1 

Dwarf 2Early 

Dwarf 2Middle 

Dwarf 2Late 

UFDl 

UFD 2 

Central 





Mvh(10«MQ) 

11.1 

7.8 

7.7 

7.6 

2.5 

1.1 

M* (10® Mq) 

620 

2200 

2700 

2800 

22 

8.5 

1st Satellite 





Mpeak (10® Mq) 

4.7 

0.70 

0.68 

0.64 

0.03 

0.02 

M* (10® Mq) 

220 

4.2 

5.4 

2.9 

(0.9) 

(0.7) 

2nd Satellite 





Mpeak (10® Mq) 

0.51 

0.04 

0.04 

0.04 

0.02 

0.007 

M* (10® Mq) 

5.0 

(0.7) 

(3.5) 

(0.6) 

(0.6) 

(0.4) 


Table 1. Halo and stellar masses of all central galaxies and their 1st and 2nd Satellites from all six simulations. For all simulations 
except Dwarf 2Barly S'lid Dwarf 2Latei the 1st and 2nd Satellites refer to the first and second most massive satellites in stellar mass. The 
2nd Satellites for Dwarf 2Barly a-nd Dwarf 2Late are selected because their counterpart within Dwarf 2Middle has M* > 3 X 10^ M©. Both 
satellites in Dwarf 1 and the 1st Satellites in all Dwarf 2 runs are “massive satellites” (see text). The stellar masses of all other satellites 
are shown in parentheses. 


mass of Mpeak ~ 5 X 10® M© our results become far less sta¬ 
ble against the stochasticity in star formation (see below). 
We therefore consider only as robust those objects that form 
within halos that have a minimnm dark matter particle num¬ 
ber Np’^ > 2 X 10® and refer to satellites that meet both this 
requirement and the stellar mass cnt as “massive satellites”. 

Dwarf 1 has two satellite galaxies that meet both of 
these requirements. Only one subhalo - the most massive 
in each of the runs of Dwarf 2 - forms a “massive satel¬ 
lite” . In Dwarf 2MiddiB there is a second satellite that meets 
the stellar mass but not the dark matter particle number 
cut. Because all three runs of Dwarf 2 use identical initial 
conditions, each subhalo in Dwarf 2Middie has a correspond¬ 
ing subhalo in the other runs that shares a majority of the 
dark matter particles. The subhalos corresponding to this 
second satellite host galaxies below the stellar mass cut in 
Dwarf 2Eariy and Dwarf 2Late. In Figure]^ stars are shown 
in this particular subhalo for all three runs of Dwarf 2. Im¬ 
portantly, although the three subhalos that host the second 
satellites of Dwarf 2 are effectively identical, they produce 
galaxies that differ by a factor of ~ 6 in stellar mass due to 
only minor changes in parameters between the runs. This 
motivates the fairly large dark matter particle number we 
require to consider a system well resolved. 

Afhaio and M* for the central galaxy (Central), most- 
massive satellite (1st Satellite), and second-most-massive 
satellite (2nd Satellite) [^in all six runs can be found in Table 
1. In the Tables, Figures and throughout the text, Mhaio is 
Mvir for centrals and Mpeak for satellites. Values in parenthe¬ 
ses indicate satellite systems that we regard as too poorly 
resolved to trust for detailed analysis. They are included 
here for completeness. The 2nd Satellites of all Dwarf 2 runs 
fall into this category, as well as all satellites of the isolated 
ultra-faints UFD 1 and UFD 2. 

In Figure the dark matter distribution is shown for 
all dark matter particles located within 65% of the virial 
radius of the halo listed in each panel. Star particles are 
only shown if they reside within any halo that currently has 


® In all cases except Dwarf 2Early and Dwarf 2Late i the 2nd Satel¬ 
lite is the second most massive satellite in stellar mass. In these 
two runs only, the 2nd Satellites are selected as the satellites that 
correspond to the high stellar mass outlier in run Dwarf 2Middle' 
Only the 1st Satellites in Dwarf 2 runs qualify as “massive satel¬ 
lites” . 


> 3 X 10® Mq in stars in the inner 3 (1.5) kpc of a central 
(satellite) galaxy or are in a corresponding halo in the other 
two runs of Dwarf 2. Many of the other subhalos also contain 
star particles (see Figure 2), but they are not shown so that 
we may focus on systems that are reasonably well resolved. 
The radial extent of the galaxies was chosen by inspecting 
their stellar mass profiles (each central galaxy also has a 
very diffuse stellar halo that we aim to avoid). The galaxy 
stellar masses are all computed using these radial limits. 


3 RESULTS 

In all six runs, the dark matter halos hosting the central 
galaxies contain significant substructure. The leftmost panel 
of Figure shows the current subhalo mass function for all 
six runs. We plot all subhalos with > 5 x 10® Mq, 

which corresponds to Vp™ > 400. Dwarf 1, with the high¬ 
est virial mass, has both the most and the most massive 
subhalos. The subhalo mass functions for the three runs of 
Dwarf 2 lie predictably on top of each other. As expected, 
the two isolated “ultra-faints” (Mvir < 3 x 10® Mq) form 
far fewer subhalos (at fixed mass) than their more massive 
counterparts. 

The satellite stellar mass functions are shown in the 
middle panel of Figurej^ They are intriguingly steep, rising 
to as many as fifty tiny satellites in the case of Dwarf 1. How¬ 
ever, the majority of these star-forming satellites form only 
a single star particle, and this potentially exciting behavior 
will need to be confirmed at much higher resolution. In or¬ 
der to make this clear, the dark green shaded region shows 
the average stellar mass of a single star particle, while the 
light green shaded region shows the mass range for all other 
objects that form with M* < 3 x 10® Mq (fewer than 16 
star particles). In all six simulations, a total of 187 subhalos 
form at least one star particle, but a mere six of them form 
galaxies with M* > 3 x 10® Mq (and only five of these also 
meet the Nf™ requirement). It is these, the most massive 
satellites, that will be the focus of our analysis (along with 
the centrals themselves). 

The M* — Mhaio relation for the five resulting “massive 
satellites” (filled triangles), their less-well-resolved counter¬ 
parts (open triangles), and all central galaxies (filled circles), 
is shown in the right panel of Figure]^ Plotted are the six 
central galaxies along with the 1st and 2nd Satellites of each 
host (detailed in Table 1) for each run. Although 18 symbols 
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Figure 2. Left: Subhalo mass function for all six runs. is the present day mass for all subhalos. The left edge of the panel 

corresponds to ~ 400 bound dark matter particles. Note: if these halos simply contained the universal fraction (/b = 0.169) of baryons, 
this would correspond to > 350 baryonic particles in these halos before they were forcefully removed by feedback. Middle: Stellar mass 
function for all satellites that form at least one star particle. The average mass of a single star particle is shown as a dark green band. 
Masses we regard as poorly-resolved (M* < 3 X 10^ Mq) are highlighted by the light green region. Right: M* — Mbaio for Centrals 
(filled circles), 1st and 2nd satellites (upward and downward triangles respectively). With surprisingly little scatter, our isolated dwarfs 
and their six well-resolved satellites (filled triangles) sit just offset from the extrapolated abundance matching relation from [Garrison- 1 
[Kimmel et al.| ) |2014| l (solid black line). Both this and the abundance-matching relationship from [Behroozi et al.| ( [^013a[ | (dotted black 
line) are shown with 0.2 dex in scatter. Open symbols correspond to poorly resolved galaxies: M* < 3 X 10'^ Mq and/or < 2 X 10®. 


are shown, we only regard the 11 filled points as reasonably 
well-resolved. All of the well-resolved systems in this work lie 
only slightly offset from the extrapolated abnndance match¬ 


ing relationship from Garrison-Kimmel et al. (20141, which 
has been constrained down to a stellar mass of ~ 10® Mq . 
This abundance matching relation, obtained by updating 


the Behroozi et al. (2013a \ relation with a faint-end slope de¬ 


rived from the measurements of the GAMA survey (Baldry 
et al. 20121, accurately reproduces the stellar mass functions 


of the Milky Way and Andromeda satellites, as well as the 
stellar mass function of the Local Field, when applied to the 
ELVIS simulations ( Garrison-Kimmel et al.|2014 |. 

We do not find the “bendin g” of the M* — Mhaio relation 
reported by Sawala et al. (20151 at Mhaio — 5 x 10® Mq . In¬ 


stead, the relation maintains a fairly steady power-law down 
to Mhaio — 10® Mq. We do see a qualitatively similar break 


to that seen by Sawala et al. (20151 but at a significantly 


lower halo mass of Mhaio — 5 x 10“ Mq. Below this point it is 
clear that we are witnessing a large amount of stochasticity, 
the onset of which corresponds to the scale of extrem ely faint 
galaxies (M* ~ 10® Mq), as opposed to the results of Sawala 
et al. (20151, where stochasticity sets in near the masses of 


classical dSphs (M* = 10® Mq). However, it is in this regime 
where our resolution is severely limited - these galaxies have 
<15 star particles. At this time we are unable to determine 
whether the apparent break and associated stochasticity in 
our relation at small mass are physical or simply artifacts of 
resolution, although it is intriguing that the break witnessed 
by Sawala et al. (20151 occurs at the mass of ~ 10 of their 


baryonic particles as well. Future work at higher resolution 
will allow us to explore this question more fully. 

The most massive satellite of Dwarf 1 represents an in¬ 


teresting statistical rarity (see Section 4.11 in that the mass 
of the subhalo hosting this satellite is just over 40% of the 
mass of its host, and the satellite itself has 1/3 of the stellar 


mass of its central galaxy. Dwarf 1 is clearly undergoing a 
major merger (see Figure and the associated large satel¬ 
lite has properties that are distinct from the other satel¬ 
lites described in this work. Most notably, its stellar mass 
is two orders of magnitude greater than the other “massive 
satellites”, and its stellar and dark matter masses are both 
larger than those of the central hosts of UFD 1 and UFD 2. 
For identification only (with no relation to the structural 
properties of the satellite), we refer when necessary to this 
most massive satellite as “dSph”, and all other satellites as 
“ultra-faint” satellites. The broader term “ultra-faints” is 
reserved for all galaxies, centrals and satellites, that have 
M* = 3 — 30 X 10® Mq in our simulations. 


3.1 Star Formation Histories 

The top panel of Figure shows the cumulative fractional 
star formation histories (SFH) for Dwarf 1 (solid line) and 
its two “massive satellites”. The SFH of its large, “dSph” 
satellite (M/^® ~ 2x 10® Mq) is shown by the dotted line. Its 
second, “ultra-faint” satellite (M/~® ~ 5 x 10® Mq) is shown 
by the dashed line. The shaded band corresponds to the 
epoch of reionization in our simulations. The bottom panel 
presents similar histories for the three Dwarf 2 runs along 
with the single most massive “ultra-faint” satellite (M/^® ~ 
3 — 5 X 10® Mq) that forms in each run (dashed lines). In 
both panels, the short vertical lines at the top mark the first 
virial crossing of the satellites. 

The massive “dSph” satellite of Dwarf 1 demonstrates a 
particularly interesting SFH. It stops forming stars shortly 
after the epoch of reionization for almost 8 Gyr (from 2 ~ 3 
to 2 ~ 0.25). However, it manages to keep some of its gas 
on hand in a state that is just below the threshold for star 
formation. Upon infall into the virial radius of Dwarf 1 at a 
lookback time of ~ 2 Gyr, followed by pericentric passage 
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~ 500 Myr later, its ISM is compressed enough to create a 
new burst of stars, all of which are concentrated at the center 
of the satellite and were formed in situ. The central galaxy 
of Dwarf 1 itself also appears to experience a related (albeit 
mild) burst. This type of bimodal star formation history was 
suggested by IRicottl] ( |2010| ) as a signature of massive reion¬ 
ization “fossils”. Our “dSph” satellite has a slightly higher 
Knax than the range given in that work, and its rebirth of 
star-formation was triggered by a merger rather than by a 
change in its halo concentration (as originally suggested by 
Ricotti|2010 1; it is nevertheless similar in its qualitative na¬ 


ture. 


In our simulations, halos with M ~ 10^° Mq appear 
to robustly maintain star formation to 2 = 0. This massive 
satellite at Mpeak — 5 x 10® Mq seems to mark a transi¬ 
tion point where reionization begins to favor uniformly an¬ 
cient star formation. The influence of its impending merger 
was enough to trigger a rebirth. The “ultra-faint” satellites 
shown in Figure]^ all form in halos smaller than this, with 
Afpeak ^ 2 X 10® Mq. They all form their stars entirely before 
2 ~ 2 and are quenched over 5 Gyr before their infall onto 
the larger hosts, indicating that their star formation is shut 
down by reionization rather than environmental processes. 
Importantly, these subsequent mergers do not trigger fresh 
star formation, as was seen in the more massive satellite. 
These small halos are either quenched during reionization 
or shortly thereafter, running out of fuel after the fresh gas 
supply was shut off by the ambient UV field. 

This tendency for our ultra-faint galaxies in Mpeak 
2 X 10® Mq halos to quench early is illustrated more clearly 
in Figure]^ where we have plotted the cumulative fractional 
star formation histories for all seven “ultra-faints” formed in 
our simulations, e.g. all “ultra-faint” satellites from Figure 
[^(dashed lines), along with UFD 1 and UFD2 (solid lines). 
We also include the unresolved 2nd Satellite of Dwarf 2Middie 
(dotted line) due to it’s high stellar mass and to illustrate 
its similar SFH. From Figure|^it is clear that, in every case, 
reionization plays a significant role in the shutting down of 
star formation in our “ultra-faints”. Not only were 100% of 
all stars in all seven of these objects formed before 2 = 2, 
four of the objects form over 90% of their stars before the 
completion of reionization at 2 = 6. 


4 DETECTING SATELLITES OF DWARFS 

4.1 How Common are Satellites of Dwarfs? 

It is beyond the scope of this paper to investigate a sta¬ 
tistically significant sample of hydrodynamic simulations of 
Mvir — 10^® Mq dwarfs in order to estimate the frequency 
with which they will host satellites above a given mass. 
However, if we assume an M* - Mpeak relation similar to 
that presented in Figure 2, we can make an estimate using 
dark-matter-only simulations. We do so using the Exploring 
the Local Volume in Simulations (ELVIS) suite of collision¬ 
less zoom-in simulations of Local Group-like environments 
( Garrison-Kimmel et al.|2014 K 

We select all isolated (dhaio > 2 for all more 

massive halos) dwarf-size dark matter halos (35 kms“^ < 
Knax < 45 kms“^) in the 12 Local Group-like pairs and the 
24 isolated ELVIS simulations, and determine the fraction of 
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Figure 3. Cumulative fractional star formation histories (SFH) 
of the Centrals (solid) and well-resolved satellites (dashed / dot¬ 
ted) in the two isolated Dwarf simulations (Mvir ~ 10^^ Mq). 
Short vertical lines along the top axes indicate infall times of the 
satellites. Top: Dwarf 1 Central (solid) along with its two satel¬ 
lites. The massive “dSph” satellite is shown by the dotted line 
and it displays a late-time burst in conjunction with its recent 
infall. The “ultra-faint” satellite (dashed) is quenched early, dur¬ 
ing reionization (shaded band), well before it is accreted. Bottom: 
Dwarf 2 centrals (solid) and the most massive (ultra-faint) satel¬ 
lites (dashed) in each run. 


those hosts that have N subhalos with Mpeak > 5 x 10® Mq. 
This is the halo mass that corresponds to M* > 3000 Mq 
according to Figure 2. ^ We similarly compute the frac- 


^ In all six of our runs, every subhalo with Mpeak > 5 x 10^ Mq 
forms a well-resolved satellite. Although we focus here on satel- 
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Figure 4. Ancient “ultra-faints”. Shown are the cumulative frac¬ 
tional star formation histories (SFH) for our isolated (solid) 
and resolved satellite (dashed) “ultra-faint” dwarfs (M^~® = 
3 — 30 X 10® Mq). The SFH of the unresolved outlier that passes 
the stellar mass cut (see Section |2.1| for details) is also shown 
(dotted line). All seven of these systems have ancient stellar pop¬ 
ulations, similar to those observed for the ultra-faint galaxies of 
the Milky Way. 


tion of dwarf halos that have a very large subhalo with 
A^peak > 4.5 X 10® Mq, Set by the most massive satellite 
of Dwarf 1 (M* > 2 x 10® Mq). 

Figurej^shows the probability for an isolated dwarf halo 
in ELVIS to have > N satellites of at least these two peak 
virial masses within a projected distance of 50 kpc (the typi¬ 
cal virial radius for our dwarfs) as a function of N. According 
to Figure isolated halos with Mvu ~ 10®® Mq will have 
one or more subhalos that could host M* > 3000 Mq satel¬ 
lites about 35% of the time. The likelihood that an isolated 
dwarf in the same mass range hosts a satellite as massive 
M* ~ 2 X 10® Mq, however, is just under 5%. 

Using both SDSS and the semi-analytic models of|Guo 


et al. (20111 applied to the MS-II simulation (Boylan- 
Kolchin et al.|2009p, [Sales et al.| (|2013[) show that the prob¬ 


ability for a central galaxy to host a satellite a given frac¬ 
tion of its stellar mass, decreases as a strong 

function of until = 10®® Mq, becoming inde¬ 
pendent of for centrals below the stellar mass of the 

Milky Way. Although the stellar masses of the central galax¬ 
ies in our simulations lie far below their stellar mass range 
(10®"'® < M*/ Mq < 10®®), this decoupling of the probability 
from for satellites of dwarf galaxies allows us to make a 
comparison with the probabilities they find. Their Figure 2 
suggests that for all central galaxies with < 10®® Mq, 
there is a ~ 40% — 50% chance that they host a satel¬ 


lites of isolated dwarfs, we note that this criteria would suggest 
the existence of ~ 100 undetected ultra-faint “stealth galaxies” 
within 400 kpc of the Milky Way. 


Figure 5. From the dark-matter only ELVIS suite of Local Group 
simulations, the fraction of isolated dwarfs in the Umax range 
35—45kms“® with N or more subhalos with Mpeak > 5x10® Mq 
(thick line) and 4.5 X 10® Mq (thin line) that are found within 
50 kpc of their host galaxy. In all of our hydrodynamic runs, all 
subhalos with Mpg^k > 5 X 10® Mq have formed a galaxy with 
Mi, > 3x 10® Mq , while the most massive satellite, “dSph”, forms 
in a subhalo with Mpg^k — 4.5 X 10® Mq. According to these 
dark-matter only simulations, ~ 35% of isolated dwarfs should 
have a “massive satellite” within their virial radii, but only < 5% 
should host a satellite as massive as “dSph”. 


life with stellar mass ~ 0.5% - similar to our 

“ultra-faints”. This is roughly consistent with the ~ 35% 
probability we find. They also find that the likelihood for a 
satellite as massive in proportion to its central as “dSph” 
(^sat/j^^cen ^ 1/3), is vanishingly small. Again, according 
to their Figure 2 and assuming that the lack of dependence 
of this probability on M/®” for dwarfs, there is only about a 
1 — 2% chance for a central to have such a massive satellite, 
which is again roughly consistent with the results shown in 
Figure 


4.2 Can they be detected? 

Our ultra-faint galaxies with M* ~ 3 — 30 x 10® Mq all 
reside within low mass halos (Mvu < 3 x 10® Mq). They 
also have extremely diffuse stellar distributions, as expected 
for the ultra-faint, low dark-matter mass, “stealth galaxies” 
discussed in BIO. Figure shows the projected (2D) half 
stellar-mass radii, R 1 / 2 , of our simulated galaxies vs their 
total stellar mass. The central galaxies are shown as colored 
circles and the satellites as triangles (the single unresolved 
satellite with M* > 3 x 10® Mq is shown as an open tri¬ 
angle). The open black circles show observed Milky Way 


dwarfs (McConnachie 20121. The black points with error 


bars are the half-mass radii and stellar masses of the eight 
recently reported ultra-faint satellite candidates from |The| 
DES Collaboration et al. (20151, and the single blue square 


with error bars is Hydra II from Martin et al. (20151. The 
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surface brightness detection limit for SDSS is shown as a 
solid black line in the Figure. It represents a constant peak 
central surface brightness for a Plummer profile 


^peak 


nRl/^ 


0.036 L 0 PC 


( 1 ) 


and corresponds to a surface brightness of fj,v = 
30 mag arcsec“^ for solar absolute magnitude M 0 ,v = 
4.83 and assuming a stellar mass-to-light ratio of M*/L = 
1(M0/L0). 

Except for UFD1, all of our simulated “ultra¬ 
faint” dwarfs have surface brightnesses fainter than 
30 mag arcsec”^, and would qualify as “stealth galaxies”. 
[^The dashed line in Figure shows a surface brightness of 
32.5 mag arcsec“^, a limit that will likely be achieved by up¬ 
coming surveys such as the Large Synoptic Survey Telescope 
(LSST). Once the full co-added LSST data is collected, satel¬ 
lites this faint should be able to detected out to ~ 1 Mpc 
over half the sky ( Tollerud et al.|2008 |. 

In the more immediate future, the Dark Energy Survey 
(DES) is already online. DES will only be able to detect 
ultra-faint objects of this kind out to about 400 kpc, and is 
only slated to cover ~ 12% of the sky over its first 5 years. 
However, DES contains within its footprint the dwarf galaxy 
Phoenix. Phoenix has a stellar mass of M* = 7.7 x 10® M© 
and is about 415 kpc away. It is fa lling into the Milky Way 
with a velocity of ~ —100 kms“® ( McConnachie|[2012 \ and 
is likely on first approach, meaning that it is a good analog to 
the isolated systems simulated in this work. This suggests 
that Phoenix may be an excellent dwarf galaxy candidate 
to host an ultra-faint satellite detectable in the immediate 
future. It also has an advantage over more distant dwarfs 
like Cetus or Aquarius in that it is close enough that low 
luminosity satellites such as the “ultra-faints” predicted in 
this work have a greater chance of being visible. 

Although, according to Eigurej^ there is only a ~ 35% 
chance that Phoenix will host a subhalo massive enough to 
form a ~ 3000 M 0 satellite, once the volume that lies be¬ 
tween us and Phoenix has been factored in, we find that 
the probability for Ending an ultrafaint in its field increases 
to anywhere between 50% — 65% according to the ELVIS 
simulations and depending on the specifics of abundance 
matching for M* ~ 3000 M 0 satellites. In general, for future 
telescopes that have a limited number of pointings available, 
targeting the ~ 50 kpc region around isolated dwarf galax¬ 
ies should prove to be a much more efficient strategy than 
pointing into blank sky, as it should increase the chances of 
observing an ultra-faint satellite by ~ 35%. 


® We have checked that the spatial radii of these galaxies are also 
well-resolved (not just their masses). Obviously, the plotted radii 
are much larger than our minimum force softening. A more de¬ 
manding criterion is the |Power et al. ]l |2003| ) radius which, for the 
dissipationless version of these runs, is > 100 pc (see |Ohorbe et alT] 
|2015| |. While this suggests that we are dynamically well resolved, 
the effect on the baryonic component within this radius is harder 
to determine. In lower resolution runs, the satellites’ half-mass 
radii do vary by about a factor of two, but in both directions 
(some are larger, some smaller). Higher resolution simulations 
would be particularly useful for solidifying the expectations for 
the sizes of these galaxies. 
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Figure 6. 2D half-stellar-mass radii, vs M* for all resolved 

galaxies in our simulations (filled circles and triangles; colors 
are the same as in Figure]^, the one unresolved satellite with 
M* > 3000 Mq (open triangle), as well as for observed Milky 
Way dwarfs < |McConnachie|2012[ open circles), the newly discov¬ 
ered DES ultra-faint dwarf candidates ( |The DES Collaboration 
et al.||2015| bla ck points with error-bars), and Hydra II ( |Ma~ 
tin et al.||2015| blue square). The solid line represents a surface 
brightness limit of 30 mag arcsec“^ while the dashed line shows 
32.5 mag arcsec“^. Both lines and the Hydra II data point as¬ 
sume a stellar mass-to-light ratio of 1. Satellites as massive as 
“dSph”, which occur around ~ 5% of isolated dwarfs (see Figure 
[^, should be currently visible, but most of the satellites lie just 
out of reach, and will visible only with future surveys such as 
LSST or potentially even DES. 


5 COMPARISON TO PREVIOUS WORK 


Sawala et al. 1 2014 1 run a series of 12 zoom-in simulations of 


Local Group-like environments at three levels of resolution, 
both with baryons and with dark matter only. They run 
each simulation with and without a cosmic UV background 
that turns on sharply at z = 11.5, and argue that the onset 
of the ionizing background radiation sets the mass scale at 
which all halos become dark. Their simulations shut down 
star formation in over 80% of halos with present-day virial 
masses less than 10® M©. Their Figure 2 does seem to show 
convergence in their highest resolution simulation run with 
reionization, but the gas particle mass in their highest reso¬ 
lution run (~ 10 ^ M 0 ) is an order of magnitude higher than 
the total stellar mass formed in most of the the satellites that 
form in our simulations, indicating that they would fail to 
detect these galaxies. Furthermore, their reionization turns 
on at a higher redshift than does ours {z = 10.65), which 
may lead to a higher halo mass cutoff for star formation. 


Shen et al. (20141 run a zoom-in simulation of a small 


group of seven dwarf galaxies with a range of virial masses 
(4 X 10® M 0 < Mvir < 4 X 10®® M 0 ). They find that galax¬ 
ies only form in the halos with present day > 10 ® M 0 , 
and of the halos that form galaxies, the two with Mvir < 
10 ®® M 0 only form stars long after the end of reionization. 
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Although they do not model H 2 cooling nor self-shielding in 
their simulations, they do approximate the effects of H 2 cool¬ 
ing and run a parallel simulation without a cosmic UV back¬ 
ground. In their reionization-free simulation, all dwarfs form 
stars, but they show that the gas in the three low mass halos 
never reaches the column density required for self shielding. 
Their baryonic particle mass is ~ 10^ Mq, the same par¬ 
ticle mass we achieve in the low resolution runs of Dwarf 2 
(see |Onorbe et S!||2015| for details). In our low resolution 
runs, each of the most massive subhalos forms an object 
with stellar mass similar to its high resolution counterpart, 
but at this resolution each satellite object consists of only 
3 — 5 star particles. While it is difficult to say conclusively 
that resolution is the main driver of their apparent halo 
mass limit for star formation, in general it has been the case 
that every time resolution is increased in simulations, the 
minimum virial mass for a halo that can form stars has ap¬ 
proached lower values | |Hoeft et al.|2006j|0norbe et al.|2015[ ). 
In addition to enabling the formation of galaxies with M* 
lower than the low resolution baryonic particle mass, higher 
resolution simulations allow for a more accurate description 
of shielding, which will also affect the minimum halo mass 
that can form stars. Furthermore, resolving the formation 
of dense substructures that collapse under self-gravity and 
become self-shielding requires resolving the Jeans/Toomre 
mass of the galaxies, which can be as low as ~ 1000 Mq 
( and corresponds to a required force/gravitational softening 
of at least < 10 — 100 pc) in baryon-poor dwarfs. 

Our simulations do not escape the challenges of imper¬ 
fect resolution, but because our particle mass is an order 
of magnitude lower than most simulations that attempt to 
form galaxies in low mass halos, we can push down the pre¬ 
dicted low mass limit for the halos that can form ultra-faint 
satellites. Rather than attempt to make a prediction for the 
lowest mass halo that can form a galaxy, we predict that 
ultra-faint dwarf galaxies can form in dark matter subha¬ 
los as low mass as Mpeak ~ 5 x 10® M©. While we cannot 
state with certainty that these galaxies are fully converged in 
our high resolution runs, the stellar masses of the satellites 
do show a better convergence between 
lution runs than do the centrals, (see 
for a discussion of the convergence of the centrals). Ad¬ 
ditionally, at these low masses, we still see sensitivity to 
spatial resolution and feedback (e.g. the second satellite of 
Dwarf 2Middie). More and higher resolution runs of dwarfs 
with Mvir ~ 10® — 10^° Mq that vary feedback prescrip¬ 
tions, reionization onset and spectrum, and particle mass 
are needed to fully probe the low mass end of the simulated 
stellar mass function and the M* — Mhaio relation. 


low and high reso- 


Onorbe et al. |[20l^ 


6 DISCUSSION AND CONCLUSIONS 


In the ACDM paradigm, all dark matter halos, from those 
around giant galaxy clusters to those hosting ultra-faint 
galaxies, should be filled with subhalos. We have used ultra- 


high resolution (m| 


255 Mq) simulations with the 


PSPH version of Gizmo |Hopkins|2014| and the Feedback in 


Realistic Environments (FIRE) prescriptions ( Hopkins et al. 
|2014[ ) to predict that subhalos of isolated dwarf halos should 
form galaxies. 

Most of the subhalos around dwarf galaxies are ex¬ 


pected to be of low-mass (Mpeak 10® Mq) and we pre¬ 
dict that some should host ultra-faint galaxies with M* < 
10'* Mq. If these tiny satellites are observed, it would pro¬ 
vide evidence that dark matter substructure persists to very 
small scales, as predicted in the standard paradigm. 

Using the dark-matter-only simulations of the ELVIS 
suite, we show that each isolated field dwarf (M* ~ 10® Mq) 
galaxy in the Local Group has about a 35% chance of hosting 
at least one satellite with M* > 3000 Mq . The extended 
~ 50 kpc regions around known field dwarfs in the Local 
Group should prove to be fruitful search areas for ultra-faint 
satellites, as each pointing towards a dwarf also contains all 
of the volume of the Milky Way dark matter halo along the 
line of sight to that region. The Phoenix dwarf galaxy in 
particular is an excellent target due to its proximity to the 
Milky Way and the high probability that it is on first infall. 

Although we consider only isolated dwarfs in this work, 
it is worth noting that some satellites of the Milky Way 
may also have their own satellites. [Deason et al.| ( |2015| ) 
show, using the ELVIS simulations, that approximately 7% 
of Milky Way satellites (but as many as 25% depending on 
the infall time and mass of the group), fell in as a part of 
LMC-size groups, and that several of the recent DES satel¬ 
lites are likely satellites of the LMC. A simple calculation 
for a rigid satellite assuming the common halo mass scale 
for Milky Way dwarf spheroidals (Mhaio = 3 x 10® Mq) 
orbiting a Milky Way-size host (Mvir = 10*® Mq) shows 
that Fornax could potentially hold on to ultra-faint satel¬ 
lites residing within its inner 17 kpc. Leo I and Leo II could 
each host ultra-faints out to ~ 30 kpc, and Sculptor out to 
10 kpc, assuming each satellite is on first infall. This sug¬ 
gests that these galaxies might also be suitable candidates 
to host “satellites of satellites”. 

We do not see a sharp cut-off or break in the M* — Mhaio 
relation, at least for Mpeak > 5 x 10® Mq. We do, how¬ 
ever, see a sharp cut-off in halos that host galaxies with uni¬ 
formly ancient stellar populations. The ultra-faint dwarfs in 
our simulations form most of their stars in the first billion 
years after the Big Bang and are subsequently deprived of 
the cold gas required for star formation due to the ionizing 
background radiation, and not by infall into a more mas¬ 
sive dark matter halo. We predict that below a critical mass 
threshold (Mhaio ~ 5 x 10® Mq, M* ~ 3 x 10'* Mq) all 
galaxies are ubiquitously ancient and, unlike more massive 
galaxies (10*” < M*/Mq < 10®) that are nearly uniformly 
star forming in the field ( Geha et al.|2012 1, both central and 
satellite ultra-faint galaxies should all be quenched. 

For satellites with M* > 10® Mq, it is likely that infall 
time is one of the primary factors in quenching their star 
formation. Massive satellites (10®'® < M*/ Mq < 10®'®) have 
been shown to have extremely long quenching timescales, 
probably due to the cutting off of the fresh gas supply after 
infall dDe Lucia et al. 2012 Wetzel et al. 2013[ [Wheeler 


et al. 20141. Slightly less massive galaxies (10® < M*/ Mq < 
10®) are likely quenched over much shorter timescales by the 
tidal or ram-pressure forces they experience upon falling in 
as satellites (Slater & Bell 2014 Fillingham et al. 20151. 
However, it is likely that the quenching of ultra-faints is 
completely independent of their infall time, and set by the 
timing of the onset of and the mass they had at reionization. 

The precise stellar masses of the “ultra-faint” dwarfs 
in our simulations are likely also sensitive to the reioniza- 
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tion redshift (with an earlier onset of reionization likely to 
quench star formation in higher mass halos) and may also 
be sensitive to star formation physics at very low metallici- 
ties. We plan to explore these dependencies in fnture work, 
but here we point out that it is possible to form ultra-faint 
satellites of regular dwarf galaxies and highlight that they 
are quenched by reionization rather than infall. 

The recent discovery of up to nine Milky Way dwarf 
satellite galaxy candidates in the Southern Sky by DES | The| 
|DES Collaboration et ^|2015| |Koposov et al.||2015| ) has 
important consequences for our understanding of star for¬ 
mation in low mass galaxies. Of particular interest is DES 
J0344.3-4331 (Eridanus II), because of its large distance 


from the Milky Way (> 330 kpc) (The DES Collaboration 
et al.|2015 l. At this distance, whether the satellite is on first 


infall or even headed out after a pericentric passage, it was 
likely accreted within the last ~ 2 Gyr. Given that the ini¬ 
tial age estimate for Eridanus II shows that it likely has an 


ancient stellar p opulation (~ 10 Gyr; The DES Collabora- 
rl this would mean that it was quenched 


tion et al. 


2015^ 


long before infall. If confirmed in follow-up observations, Eri¬ 
danus II would be the first known ultra-faint galaxy shown 
to be quenched in the field, supporting our findings that, 
below a critical mass scale ~ 5 x 10® M©, all galaxies host 
ancient stellar populations quenched by reionization-related 
feedback, and not by environmental processes. This would 
be a clear example of a reionization “fossil”, as first discussed 
by Ricotti & Gnedin (20051. 
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